Analogue black hole in magnetohydrodynamics 
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We consider an irrotational plasma fluid evolving under the effect of a background magnetic field. 
The magnetohydrodynamic formalism is used to describe the electromagnetic waves and the dynam- 
ics is described by a scalar field that follows a second order differential equation. This equation can 
also be recovered as the wave equation associated to a field in a curved space-time. Through this 
analogy we recreate a sonic horizon, equivalent to those found in perfect fluid theories. However, 
in this case, the magnetic field creates a pressure in the plasma which contributes to the magne- 
toacoustic speed that builds the horizon. This effect enhances the temperature produced by the 
Hawking radiation expected from this analogue black hole, and eventually, making its experimental 
detection worth to consider. 
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I. INTRODUCTION 

Thirty years after Unruh's seminal paper [l[ , the study 
of the field of artificial black holes is attracting the inter- 
est of a large number of researchers. The result found by 
Unruh emerges from the similarity between fluid theory 
and general relativity at the kinematical level. Those 
systems whose dynamics is controlled by a second or- 
der differential equations display the same features than 
a dynamical field evolving on a curved space-time. In 
these cases, an effective metric can be extracted from the 
equation of motion of the fluid, and this metric display 
similar properties to those found in the neighborhood of 
a black hole event horizon. 

A black hole exhibits what is called the Hawking ra- 
diation which is associated to the event horizon with a 
temperature proportional to the surface gravity of the 
black hole @- However, the detection of this radiation is 
a real challenge. In this context, if the analog black holes 
may reproduce all the features of the Hawking radiation, 
it offers the opportunity to search for the quantum insta- 
bility proposed by Hawking. For this reason, there are 
several ongoing efforts searching for this radiation. The 
experimental verification and the possibility of studying 
their properties within this new context, would show the 
robustness of the Hawking effect. 

Some of the settings for analog black holes include, per- 
fect fluids, where Unruh originally found a horizon gener- 
ated by the relative changes between the acoustic and the 
flow speeds P10|. This effect has also been studied within 
the context of gravity waves Q and dielectric mediums 
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Q , looking for systems where a higher Hawking temper- 
ature could be measured. The search has been concen- 
trated on the Bose-Einstein condensates [fUl], where it 
has been claimed a measured evidence of the Hawking ra- 
diation (Io| . In this context, they estimated a spectrum 
of phonons at ~ 7 • 10~ 8 K @. The use of waveguides 
have also been proposed to detect the Hawking radiation 
since the electromagnetic radiation can be controlled and 
detected much easier than sound [ll[ . Using this mecha- 
nism, the experimental formation of an artificial horizon 
have been achieved in optical fibers [l2j . 

Recently, an ultrashort laser pulsed experiment has 
been performed to verify the Hawking radiation (l3| . The 
results are encouraging even though they have rose some 
controversy concerning the thermal spectrum of the radi- 
ation measured flij . Another experiment using an open 
channel fluid with an obstacle interposed [f5[ has allowed 
to mimic a white hole, verifying the thermal nature of the 
emitted radiation, as underlined by the authors. 

Here we address the occurrence of analogue black holes 
in plasmas. We compare the Hawking effect associated 
to sound waves in this context with the known ones using 
neutral fluids. Magnetized plasmas offer a new kind of 
contribution to the analogue black holes, since they are 
the simplest system where a fluid of charged particles 
interact with a magnetic field. 

Concretely, we study a magnetized plasma using the 
ideal magnetohydrodynamics (MHD) theory. In this ap- 
proach, a plasma, composed by electrons and ions, is 
described as a single fluid. The equations for MHD are 
constructed as a linear combination of the fluid variables 
of the two kind of charged particles [l6j . 

For a quasincutral plasma composed by electrons and 
ions, a one-fluid mass density can be defined as p = 
n (M + m), where n is the density number, equal for elec- 
tron and ion fluid, and M and m are the ion and electron 
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masses respectively (M ^> m). Similarly, a single fluid 
velocity is defined as v = (Mv, + mv e )/(M + m), where 
Vj and v e are the velocities of the ion and electron fluid, 
respectively. 

When the Hall current effect is neglected and a null 
resistivity plasma is considered, it is possible to find a 
simple set of MHD equations, called ideal MHD. In this 
case, the Ohm law is written in the simple form E + v x 
B = 0, and the displacement current d~E/dt is neglected 
from the Ampere's law. 

The equations that enter in this description for the 
ideal MHD theory are, the continuity equation 



dp 

at 



V ■ (pv) = . 



(1) 



the equations of motion for a single fluid 



p(| +v .v)v = -v(p + i|Bp) + i-(B.V)B, 

_ (2) 

where p is the pressure of the single fluid [17[, and 
Maxwells equations are 

— = (B-V)v-B(V-v)-(vV)B. (3) 

Eqs. @ and are obtained using the current density 
J = (c/47r)V x B, from Ampere's law, and the Ohm law 
respectively. 

We can use the ideal MHD equations to find the prop- 
agation of waves in any configuration of electromagnetic 
and fluid fields. 



II. MAGNETOHYDRO DYNAMICAL BLACK 
HOLE 

Since we follow a perturbation scheme, we define a 
subindex for the non-perturbed background fields, and 
a subindex 1 for the first order perturbations. So, p = 
Po + Pi, v = v +vi, B = B +Bi andp = p +pi- The 
background magnetic field Bo and the background den- 
sity po remain constants in this context, and consequently 
the background becomes incompressible, V ■ v = 0. 

In general, the magnetic field introduces a vorticity 
into the plasma. However, the longitudinal modes ac- 
cept an irrotational solution for the plasma velocity. The 
magnetoacoustic mode is one of these longitudinal modes 
where the wave is longitudinal in the perturbation of the 
fluid velocity, although the perturbation of the magnetic 
field is transverse. Considering perturbations of the form 
cxp(ik-x — iut) in the perturbed quantities, where k and 
uj are the wavenumber and its frequency respectively, the 
propagation is transversal to the perturbed field Bi, i.e. 
Bi • k = 0. This condition is relevant for the propagation 
of magnetoacoustic modes. 

We also impose the following constrains: the perturba- 
tion Bi is chosen to be parallel to the background mag- 
netic field Bp x Bi = 0, and then Bp- k = 0. This implies 



that (Bo • V)Bi = 0. Besides, we choose Bo • vo = for 
the longitudinal mode. Given these assumptions, we pro- 
ceed to write Eqs. ((J)- ([3]) keeping only first order terms 



d_ 

at 



v • V pi + p V • vi = . 



(4) 



<9 V 1 t-, / ^ 

Po-^-+poV(v -vi 

at 



Vpi - — V (B ■ B0 , (5) 
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— +v - VJ Bi = (Bi • V)v +(B - V)vi-B (V-v0 

(6) 

where the perturbed velocity is longitudinal. The Ohm's 
law is satisfied to first order: Ei + vi x Bo + vo x Bi = 0. 

From Eq.© we have that Vxvi =0, the perturbed 
velocity aligns itself in the direction of the wave propa- 
gation k and, as a consequence, the plasma fluid remains 
irrotational. This is the main step behind our approach. 
When a fluid is irrotational, its velocity can be described 
by a single scalar field potential ip\. The first order per- 
turbation velocity is then, determined by this potential: 
Vl = Wi- 

For the plasma to remain irrotational as it evolves in 

time, the pressure associated to the fluid should depend 

only on the density [l8|. In our case, we have chosen 

an isothermal plasma, such that V pi = v 2 Vpi , where 
l li 

v s = (7P0/P0) is the adiabatic sound speed. When 
ions and electrons are involved in a fluid, the sound 
speed is given by v 2 = k B ( 7» T i + le T e )/(M + m) f« 
ks {liTi + jeT^/M, where Ti and T e are the ion and 
electron fluid temperatures respectively and ks is the 
Boltzmann constant [Tgj | . In this case, we use 7 = 5/3. 

With this choice for the magnetic fields and with the 
condition of an irrotational plasma, we have that Bo • 
vi = . Then, from equations (Q| and Eqs.©, we get 



d_ 

dt 



vo 



V (Bo-Bi -65 2 ) =0. 



(7) 



where £i = pi/po and Bq = B • B . Therefore, using 
(J?]) and neglecting the integration constant, the Eqs. (0| 
and ([5]) become 



d 



v a + v> x = o , 



(8) 



— + v • V ) i> x = -c 2 ^i , (9) 



where c is the magnetoacoustic speed [2Cj defined as 

,V2 



and 



Va 



B 



V 47r /°o ' 



(10) 



(11) 
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is the Alfven velocity. 

The equations © and © for ^ and ipi show the mu- 
tual influence -as expected- from the dynamics of the 
density and the velocity perturbations for the irrotational 
plasma fluid. The magnetic field is proportional to the 
Alfven velocity and, through this velocity determines the 
magnctoacoustic speed, see (|10p and (fTTj) . The lineariza- 
tion of both equations gives the propagation of the mag- 
netoacoustic modes in the MHD formalism [20|. This 
mode propagates perpendicularly to the magnetic field 
compressing and releasing both the lines of force and the 
conducting fluid. 

Eqs. ([8]) and can be combined into a single second- 
order differential equation for the scalar field ipi 



at \P^ + p^o-y^ 



'dip! 



p Q v ■ V 

-p c 2 V 2 V'i = , 



i v • V^i 

ot 



black hole as first proved by Unruh Black holes are 
characterized by an event horizon, a mathematical null 
surface that splits the trajectory of light rays depending 
if they travel inside or outside this surface. This is a 
kinematical consequence of the event horizon existence. 
On the other hand, the entropy as proportional to the 
area of a black hole is a dynamical consequence of Ein- 
stein equations [2l[ . In the analogue black hole, the event 
horizon is a kinematic effect and is characterized for the 
velocity of the fluid reaching the sonic velocity. It is ex- 
plained in more detail with a drawing in the next section. 
The temperature arises in the analog black hole from the 
quantizing of the scalar field characterizing the velocity 
perturbation. This is the main result from the first paper 
by Unruh The procedure is generic and it is possible 
to follow step by step for this case. We return to this 
subject in the following section. 



(12) 



where the coefficients are given by the background quan- 
tities of the plasma. This equation contains the same 
information about the propagation of the magnctoacous- 
tic waves than the two previous equations. 

Eq. (|T2|) can be re-arranged to appear as the Klein- 
Gordon equation for a masslcss scalar field in a curved 
space-time: <9 M (yj—gg^ d„ipi) = 0- To achieve 

this, we need to introduce the effective Unruh metric [l[ 
(with the Greek indices running from 0-3) 



1 



-1 



(13) 



where v J is the j component of v (the Latin indices run 
from 1-3), and g = [detg^]^ 1 = — c 2 . The inverse of 
this metric is 



Po vfi-c 2 



(14) 



where t> 2 = vo ■ vq. 

The metric (fbfj) can be interpreted as an effective 
curved space-time where these waves propagate. The in- 
terval extracted from this geometry is 

ds 2 = g^dx^dx" 

= — (N) - c 2 ] dt2 - 2dtv o ■ dx + dx • dx) (15) 

Assuming that the background flow is spherical sym- 
metric and stationary, and defining the new time dr = 
dt — vodr/ (vq — c 2 ) [l(, we can write the interval (|T5|) as 



ds 2 = ^ ( (v 2 - c 2 ) dr 2 



r 2 dn 2 



(16) 



where dfl 2 = dO 2 + sin 2 ddtp 2 . In general, vo depends on 
the radius. 

The metric (|16p assuming a velocity vq ~ — c + a (r — 
R) + ... displays properties similar to a Schwarzschild 



III. HAWKING RADIATION IN PLASMAS 

The possibility of the existence, in the sense of mea- 
surable, of a thermal radiation coming from an acoustic 
horizon caused by a supersonic flow, will set the thermal 
radiation generated on the event horizon of a black hole 
in a different perspective. For this reason we examine 
the Hawking temperature obtained for a MHD plasma 
within this context. 

In Unruh's approach [l[ there exists a scalar field ip, as 
a perturbation around some background solution ipo re- 
lated to the fluid velocity potential. Expanding this field 
tjj in terms of creation and annihilation operators and af- 
ter some canonical manipulation he concluded that for 
this fluid with a sonic horizon, there exists a thermal 
spectrum of sound waves radiated. The expression ob- 
tained for the temperature of the analogue black hole can 
be compared term by term with the Hawking expression 
and given a corresponding physical interpretation. 

Following the same steps applied to the MHD fluid, 
it occurs that near the acoustic horizon, the analogue 
black hole will emit magnetoacoustic waves with a ther- 
mal spectrum with a temperature given by 



T H = 



hgn 
2tt ks c 



(17) 



where gu is the surface gravity on the acoustic horizon 
[l8| given by the following expression 



-— (C — Vq) 

an 



(18) 



and d/dn represents the derivative of the plasma fluid 
velocity taken normal to the cross section of the horizon. 
In our case, c > v s always due to the expression of the 
Alfven velocity. Thus, the background magnetic field, 
within this context, increases the value of the Hawking 
temperature compared with the effect of a neutral fluid. 
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We propose the following setting to achieve a magne- 
toacoustic wave propagating through a background mag- 
netic field and its thermal radiation associated, as de- 
scribed in Sec. |TIJ Let's consider a plasma traveling 
in a cylindrical configuration under the influence of a 
background magnetic field in the transversal direction 
Bo = Bqz (see FigfT]). These waves will propagate in 
the direction of increasing 9. We need to introduce an 
approximation at this point, the radius of this cylindri- 
cal configuration must be large compared with the size 
of the nozzle. In this way when we compute the relevant 
physics in this region assuming a straight motion of the 
plasma. 

Following Ref. [l8[ , we set a Laval nozzle in some region 
of this cylindrical device. Within this device, when the 
cross section shrinks, reaching its minimal section, the 
plasma flowing through it, will go supersonic. For this 
condition to occur, as mentioned above, the radius of the 
cylindrical device must be larger than the length charac- 
terizing the nozzle. In other words, the Larmor radius of 
the plasma constituents induced by the background and 
perturbed magnetic fields must be much larger than the 
lenght of the Laval nozzle. Within this approximation 
the plasma can be considered as a longitudinal flow in 
that region. In FigQ] we show a diagram of the Laval 
nozzle and the background magnetic field. 

If the space-time dependence of this wave is exp(ik • 
r — iust), where k = k9, then from Eqs. (JSJ and ©, we 
can obtain the dispersion relation for the magnetoacous- 
tic waves v g = »o ± c [2Q |. where v g = dui/dk is the 
group velocity and the drift velocity of the plasma lies 
mainly on the longitudinal direction Vo = vq9. It is easy 
to see that if the fluid is subsonic, then v g < for the 
reflected waves in the fluid and v g > for those trans- 
mitted, therefore the waves travel forward or backward 
along the 9 direction under this condition. At the hori- 
zon, the center of the nozzle, when the velocity becomes 
supersonic vq = c, occurs that v g > for all waves, and 
every wave crosses the horizon without return. This noz- 
zle behaves as a black hole event horizon. At the other 
side of the nozzle, those waves become supersonic with a 
group velocity v g > and v g > c always, and they keep 
traveling forward (to the right). 

A time inverse situation occurs in the following nozzle, 
as the velocity decreases, reaching vq = c at the noz- 
zle and there, the velocities are equal or larger than the 
sound speed and all the fluid moves to the right. It is the 
equivalent of a white hole. 

The Laval nozzle is the key to generate a supersonic 
plasma fluid. It plays the role of the acoustic horizon. If 
the radius of the narrowest part of the nozzle is R, we 
can approximate 



1 d(c 
c 



dn 



1 

R 



(19) 



This approximation is valid because the magnetoacous- 
tic speed is constant. Then, the surface gravity will be 




v 0=c 
v g>0 



FIG. 1: Laval nozzle and background magnetic field. First, 
the wave is subsonic with v g < or v g > 0. The narrowest 
part of the nozzle represents the horizon with vo = c and v g > 
(the black hole, BH). Then, the waves become supersonic 
with v g > only. Finally, they reach the other horizon (the 
white hole, WH) to return to be subsonic. 



9h ~ c 2 /R. This is essentially the same result obtained 
for the surface gravity gn ~ c 2 / \J2 Ah considering the 
cross section area Ah ~ R 2 of the Laval nozzle @. 

Thus, applying the same formalism used in previous 
works [lj, [6j , the corresponding Hawking temperature is 



Ti 



he 



H 



2nk B R 



(20) 



This is the main result of this work. Also, when we 
compare this temperature with previous cases, it ap- 
pears that we always can achieve a larger temperature 
compared with a neutral fluid because c > v s . An in- 
teresting case occurs when the Alfven velocity is much 
larger than the sound speed. The temperature of a ten- 
uous laboratory plasma is around 10 K, and therefore 
the sound speed will be of the order of 10 6 cm/s. On 
the other hand, for typical values of the number density 
of no = 10 12 cm -3 and a magnetic field of £>o = 10 4 
gauss, the Alfven velocity is of the order of 10 9 cm/s, 
and hence » s < b^. In this case, the magnetic pressure 
becomes larger than the fluid pressure B 2 /4ir 3> po, and 
then c ~ va- The Hawking temperature now will de- 
pend mainly on the background magnetic field through 
the Alfven velocity. Recalling that m/M ^s> 1, then 
po = M no- The Hawking temperature becomes 



T H « 2.66 



1/2 /l mm \ (\ cm" 3 \ 1/2 / B a \ 

V R J V n o J \g auss / 

(21) 

where m p is the proton mass. 

It is possible to imagine configurations where the mix 
of a magnetic field and an appropriate number density 
can produce a Hawking temperature larger than that as- 
sociated to a neutral fluid. Keeping in mind that the size 
of the nozzle must follows the geometric limit imposed 
for our approximation. 

With the following estimations: the ion mass as the 
proton mass and the radius of the nozzle considered as 
1mm, then for the number density equal to uq = 10 12 
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cm" 3 and with a magnetic field Bq = 10 4 gauss, we ob- 
tain Tjj ps 0.027 K. However, the Hawking temperature 
can grow higher with other values for larger magnetic 
field and for hotter plasmas. 



IV. CONCLUSIONS 

Here we have studied and concluded that a plasma dy- 
namics under the magnetohydrodynamics approximation 
mimics closely certain aspects of scalar fields in curved 
space-time just as a normal fluid docs. The relative ve- 
locity associated to the magnetoacoustic modes and that 
of the flow itself, generates an analogue horizon equiva- 
lent to the event horizon of a Schwarzschild black hole. 
The similarity between these two physical systems is ro- 
bust. The Hawking temperature generated by the event 
horizon it is replicated in this plasma evolving on a back- 
ground magnetic field chosen adequately. The Alfven 
velocity increases the magnetoacoustic speed as shown 
in (|10j) . and since this velocity depends on the external 
magnetic induction Bq, it is feasible to raise it. The den- 
sity of the plasma is another parameter that affects the 
Hawking temperature in this plasma. With an increasing 
density, the hawking temperature decreases. 

Due to the the contribution of the magnetic field then, 
the equivalent of the Hawking radiation in a plasma be- 



comes larger than that in neutral fluids, and it can even 
grow higher for larger magnetic fields and temperatures. 
However, the magnetic induction cannot grow arbitrarily 
large since there are geometric constrains as pointed in 
the text. Also the plasma requires temperature as large 
as 10 4 to 10 5 degrees. So the ratio between the temper- 
atures of the the thermal emission and the background 
decreases moderately, with a factor of 10 6 or 10 7 . Nev- 
ertheless in neutral fluids the ratio is larger, around 10 9 . 

In this work we have focused ourselves only in the gen- 
eral setting of the problem. We have not given a specifics 
of such a system. However, the most recent published 
results (considering other fluids) give a feeling that this 
effect is robust and remains detectable directly or indi- 
rectly even though the external conditions are not the 
ideal, see for instance S. Weinfurtner et al. 15[ where the 



thermal nature of the emission appears to be confirmed 
and also in the work of Belgiorno et al. [l3[, where the 
radiation detected, according to the authors, can only 
be attributed to an equivalent of the Hawking radiation. 
There are more experiments proposed, for instance in 
Ref. [HI, H3. All these results make feasible that the 
effect displayed here can be possible measured. 

All these elements make this problem more interest- 
ing, worth to be studied further. Also, the simplicity 
and ingenuity displayed in the experiments, contributes 
strongly in this direction. 
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